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rH* ■ ABSTRACT 

a.: 

The prototypical accretion-powered millisecond pulsar SAX J1808.4-3658 was observed simultaneously with Chandra-LETGS and 
RXTE-PCA near the peak of a transient outburst in November 201 1. A single thermonuclear (type-I) burst was detected, the brightest 
yet observed by Chandra from any source, and the second-brightest observed by RXTE. We found no evidence for discrete spectral 
features during the burst; absorption edges have been predicted to be present in such bursts, but may require a greater degree of 
photospheric expansion than the rather moderate expansion seen in this event (a factor of a few). These observations provide a unique 
data set to study an X-ray burst over a broad bandpass and at high spectral resolution (A/AA = 200-400). We find a significant 
' excess of photons at high and low energies compared to the standard black body spectrum. This excess is well described by a 20-fold 

increase of the persistent flux during the burst. We speculate that this results from burst photons being scattered in the accretion disk 
^v^j , corona. These and other recent observations of X-ray bursts point out the need for detailed theoretical modeling of the radiative and 

■ hydrodynamical interaction between thermonuclear X-ray bursts and accretion disks. 

Key words. Accretion, accretion disks - X-rays: binaries - X-rays: bursts - stars: neutron - X-rays: individual (SAX J1808.4-3658) 

t-H . 

< ^ > 1. Introduction Paizis et al. 2012). Hardly any X-ray bursts with photospheric 

expansion and low A^h were detected with sub 2 keV coverage 
. . ■ The vast majority of thermonuclear X-ray bursts from neutron so far, perhaps 6 with low flux (Galloway et al. 2010). 
> ; stars (the so-called type I X-ray bursts, see e.g. Lewin et al. ln m i s paper we report an observation of a bright thermonu- 
1993; Strohmayer & Bildsten 2006) have been measured with c i ear burst with photospheric expansion from the low-mass X- 
instruments that are sensitive at photon energies above 2 keV. ray binary (LMXB) SAX J1808.4-3658, with the AXAF CCD 
; The temperature of type-I X-ray bursts is, however, 2.5 keV or i mag ing Spectrometer (ACIS) in combination with the Low- 
lower so that usually a substantial part of the spectrum is missed. Energy Transmission Grating Spectrometer (LETGS) on board 
Therefore, interesting physics may be missed, particularly dur- Chandra. The event was detected simultaneously with the RXTE 
ing the cooler parts of X-ray bursts in the tail and during phases Proportional Counter Array (PCA). Combined this presents a 
of photospheric expansion. Not only that: many of the discrete unique data set with, as far as we know, the largest number of 
atomic spectral features from abundant elements in the cosmos photons ever detected from a single X-ray burst in the 0.5-2 keV 
(oxygen, neon and iron) occur below 2 keV. Sub-2 keV mea- band or with a grating spectrometer. 

surements of X-ray bursts at high spectral resolution and with SAX ji 808 .4-3658 is the first discovered accretion- 

photosphenc expansion are, therefore, interesting to pursue. p0W ered millisecond X-ray pulsar in a LMXB (in 't Zand et al. 

Since the launch of Chandra and XMM-Newton in 1999 and 1998; Wijnands & van der Klis 1998; in 't Zand et al. 2001; 
Swift in 2004, sub 2 keV coverage is readily available. However, Chakrabarty et al. 2003), of 14 cases known thus far (e.g., 
these are narrow field instruments where measurements are usu- Patruno & Watts 2012). It is also a relatively nearby represen- 
ally made through a dedicated program, in contrast to the wide- tative, at 3.5 ±0.1 kpc (Galloway & Cumming 2006), and with 
field instruments on board SAS-B, BeppoSAX, INTEGRAL, a relatively low column density of interstellar material (A^h = 
Swift and Fermi that lack sub 2 keV coverage. We estimate (1 - 2) x 10 21 cm" 2 ; e.g. Wang et al. 2001; Papitto et al. 2009; 
that about 250 X-ray bursts have been detected below 2 keV Cackett et al. 2009). Almost all its X-ray bursts exhibit photo- 
with these three missions, compared to about 10 4 with the wide spheric expansion (Galloway et al. 2008) and the bursts appear 
field instruments (e.g., Cottam et al. 2002; Jonker et al. 2003; to be the result of a flash of a pure helium layer that is pro- 
Thompson et al. 2005; Boirin et al. 2007; Kong et al. 2007; duced by stable hydrogen burning. These circumstances make 
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Fig.l. Light curves of the 2011 outburst of SAX J1808.4-3658: 
from 2-10 keV RXTE/PCA data (green dash-dotted line, normalized 
to 1580.5 c s^'PCLT 1 equivalent to 1.0 Crab unit; see also Patruno 
et al. 2012), 2-10 keV MAX1/GSC data (red dashed line; normalized 
to 3.2 c s~'cirr 2 equivalent to 1.0 Crab unit) and 15-50 keV SwiftfBAT 
data (black solid line; normalized to 0.22 c s~'citT 2 equivalent to 0.18 
Crab units). The outburst started on MJD 55865 (Oct. 31st) and lasted 
for 21 d, until MJD 55886. The horizontal solid line indicates the time 
frame of the Chandra observation and the vertical dotted line through it 
the burst that was detected with Chandra and RXTE. The vertical bars 
at mid left indicate the typical l<x uncertainty per data set (PC A flux 
errors are typically smaller than the line thickness). 



SAX J1808.4-3658 a particularly suitable source to study X- 
ray bursts in great detail and in a large bandpass. Therefore, we 
observed SAX J1808.4-3658 with Chandra and RXTE during 
its November 201 1 outburst. 

The main goal of our measurement was to search for absorp- 
tion edges. We were motivated by a prediction of Weinberg et al. 
(2006) that the radiative wind of radius expansion bursts eject 
ashes of nuclear burning whose spectral signature may be de- 
tectable with current high-resolution X-ray telescopes. This may 
be the best path to finding discrete spectral features from neu- 
tron star (NS) surfaces, whose detection could constrain ther- 
monuclear reactions in the NS ocean and the NS compactness 
through general relativistic gravitational redshift determinations. 
Near-Eddington X-ray bursts reveal themselves as having black 
body like spectra with phases of expanded emission areas. There 
are two kinds: those with moderate expansion (factor ~ 10 in- 
crease in emission area), encompassing « 20% of all bursts, and 
those with superexpansion (factor ^ 10 4 ), encompassing perhaps 
~ 0.1% of all bursts (e.g. van Paradijs et al. 1990; Galloway et al. 
2008; in 't Zand & Weinberg 2010). The difference in expansion 
is likely determined by whether an optically thick shell is being 
expelled which on its turn is determined by how much the radia- 
tive flux from the nuclear burning goes over the Eddington limit 
in the flash layer. Also, the larger this super-Eddington factor 
is, the larger the upward extent of the convection and the prob- 
ability for absorption edges will be. The burst we discuss here 
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Fig. 2. Time histories of total C/iandra-observed photon rate, measured 
over all active ACIS-S CCDs, at 100 s resolution (top panel), the color 
at 1000 s resolution (middle panel) and the RXTE PCU2 rate at 100 s 
resolution (bottom panel). The burst peak has been cut in all panels. No 
background was subtracted. 



is a bright burst from a relatively nearby NS and, ergo, of good 
statistical quality. However, it has only a moderate photospheric 
expansion. 

We introduce the observations in Sect. 2, including moni- 
toring observations with the Monitoring All-sky X-ray Imager 
(MAXI) on board the International Space Station (ISS) and the 
Burst Alert Telescope (BAT) on board Swift, and provide the de- 
tails of the data reduction in Sect. 3. Sect. 4 discusses the spectral 
analysis and Sect. 5 briefly the timing behavior of the detected 
burst. Finally, in Sect. 6, we place the results in the context of 
the theory of burst spectra. 

2. Observations 

On November 4, 2011, a new outburst of SAX J1808.4-3658 
was reported from Swift-BAT observations beginning October 3 1 
(Markwardt et al. 2011). Fig. 1 shows the outburst light curves 
from MAXI, RXTE-PCA and Swift-BAT data. MAXI (Matsuoka 
et al. 2009) is deployed on the International Space Station (ISS) 
and, in combination with the Gas Scintillator Counter (GSC), 
is active between 2 and 30 keV. The measurements consist of 
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Fig. 3. Time history of RXTE detected photon rate around the burst, at 8 
s resolution. Photons detected with PCUs 1 and 2 have been combined. 
The dotted horizontal line shows the fit with a constant line to the data 
prior to the burst (with a best-fit value of 418 s -1 ). The burst lasts 550 + 
20 s. The inset zooms in on a 25 s stretch around the peak, at 1/8 s 
resolution. 



1 min exposures every 100-min ISS orbit and are made public 
within a few hours of being taken 1 . BAT (Barthelmy et al. 2005) 
is a gamma-ray burst monitor camera active on Swift between 15 
and 100 keV. The 15-50 keV flux of some 500 X-ray sources are 
monitored several times per day and made available publicly in 
near real time 2 . The PCA data (see below) were obtained under 
ObsID 96027-01 for a total exposure time of 143.6 ks (see also 
Patruno et al. 2012, and below). 

Following the alert, we requested and were granted a 
Chandra Target of Opportunity Observation (ObsID 13718). 
The exposure time is 100.13 ks between Nov. 8 00:46 UT and 
Nov. 9 04:35 UT (see Fig. 2) which is 8 days after the presumed 
outburst onset. The observatory was operated using the Low- 
Energy Transmission Grating Spectrometer (LETGS; Brinkman 
et al. 2000), with ACIS-S (Garmire et al. 2003) at the focal plane 
and the Low-Energy Transmission Grating in the beam. This 
combination provides a marginally higher effective area than 
the HETGS (see e.g. Proposers' Observatory Guide Fig. 1.4) in 
most of the interesting bandpass of 0.5-6 keV at the expense of 
some spectral resolution. The resolution is, however, sufficient to 
search for absorption edges (A/AA = 200-400). A disadvantage 
of using ACIS instead of HRC is the loss of effective area be- 
tween 0.08 and 0.5 keV, but this bandpass is of no interest to this 
observation because the source is significantly absorbed there. 
To avoid pile up, ACIS was operated in the continuous clock- 
ing (CC) mode. In this mode, position information is lost in the 
cross dispersion direction while the time resolution is improved 
from 3.2 s to 2.85 ms. The disadvantage is that the background 
cannot accurately be measured independently (it is possible to 
discriminate with some sensitivity the background from source 
in pulse height versus dispersion space), but this is not detrimen- 
tal because we are primarily interested in X-ray bursts which can 



1 http://maxi.riken.jp/top 

2 http : //swift . gsf c . nasa . gov/docs/ swift/results/ 
transients 



easily be separated in the time domain, without the need for the 
spatial resolution. Four ACIS CCDs were employed: S 1 through 
S4. The offset of the Science Instrument Module (the so called 
'SIM offset') was -8 mm to position the dispersed spectrum as 
close as possible to the CCD edges to optimize the CCD spectral 
resolution. The Y offset was + 1 .6 arcmin so that the zeroth order 
would be detected by S2. 

Simultaneous measurements were requested and granted 
with RXTE under Proposal Number 96027. Twenty-four ks of si- 
multaneous coverage was obtained (see Fig. 2). The Proportional 
Counter Array (PCA; Jahoda et al. 2006) comprises 5 non- 
imaging proportional counter units (PCUs) active between 2 and 
60 keV with a combined effective area of 6000 cm 2 at 6 keV. 
The spectral resolution is 18% full width at half maximum at 
6 keV and the time resolution of the data products for SAX 
J1808.4-3658 is typically 125 fis. It is seldom that all PCUs are 
active at the same time, particularly at late times in the mission. 
During our observations, PCA operated with 2 active propor- 
tional counter units (PCUs 1 and 2, counted from 0). All PCUs 
have collimators delimiting the field of view to 2°x2° (full width 
to zero response). For our analysis, we used standard- 1 and 2 
data (time resolution 0. 125 and 16 s, energy resolution 1 and 128 
channels, respectively) and event mode data E_125us_64M_Q_ls 
with 125 //s time resolution and 64 energy channels between 2 
and 60 keV. All data allow selection of PCUs. 

Fig. 2 shows the Chandra and PCA-measured (only PCU2) 
light curves during the time of the Chandra observation. These 
curves show the detection of one burst in both instruments. 
This is the only X-ray burst detected from SAX J1808.4-3658 
in the whole outburst (c.f., Patruno et al. 2012). RXTE and 
Chandra times were synchronized by matching the burst on- 
set as seen in the overlapping bandpass 3-8 keV. We find 
that the burst starts at MJD 55873.91320 (barycentered), or at 
RXTE time 3 563407175.90s and Chandra time 4 437176836.46s 
(not barycentered). There is simultaneous Chandra/RXTE cov- 
erage between -2560 and +812 s with respect to the burst on- 
set time (see Fig. 3). The total number of burst photons de- 
tected is 38,800±180 with Chandra (measured over the com- 
plete ACIS-S detector and after subtraction of a pre-burst level) 
and 562,800±400 with RXTE. The peak intensities, measured at 
0.125 s time resolution, are 2,538 c s _1 for Chandra (all ACIS 
events, not corrected for non-burst events) and 33,669 c s _1 for 
RXTE. The Chandra net peak rate for orders - 1 /+ 1 is 1 ,930 c s~ 1 . 
The Chandra event data buffer of 128,000 events (in the graded 
telemetry mode) was not overfilled. The RXTE data buffer was 
overfilled. Starting at 1.15 s after burst onset, 0.45 s chunks of 
data are missing in the event-mode data every 1.0 s until 10.6 
s. After that 3 smaller chunks of 0.45 s are missing every 1.0 s 
until 13.6 s after which the data are completely recovered. The 
standard- 1 and 2 data do not have gaps. 

The non-burst emission varied significantly during the 
100 ks of the Chandra observation. The intensity decreased from 
128 to 110 c s _1 , or 14% although this is not corrected for the 
background (which is expected to be a few c s~'). The decrease 
is not smooth. The Chandra color history, the color being de- 
fined as the intensity of all ACIS events above 1 keV divided by 
that below, also shows a gradual decrease by 5%. 



3 Mission Elapsed Time, defined as seconds since 1994.0 UTC 

4 Mission Elapsed Time, defined as seconds since 1998.0 TT 
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Table 1. Spectral parameters of fit to non-burst or persistent 
Chandra/RXTE spectrum in 2.5 ks prior to burst. 



Parameter 



Fitted value 



disk black body kT m 

black body norm. (&J D 2 l0kpc ) 

power law photon index 

power law flux at 1 keV 

gauss centroid 

gauss sigma 

gauss normalization 

xl/v 

0.5-20 keV flux 
0.5-60 keV abs. flux 
0.5-60 keV unabs. flux 



(0.! 



: 0.04) x 10 21 cnT 2 



0.59 1 ±0.006 keV 
(6.5 ± 0.3) x 10 2 
1.892±0.013 

0.399±0.014 ph s-'cm^keV- 1 
5.12+0.09 keV 
1.66+0.06 keV 
0.0156±0.0012 phot s^'cnr 2 
0.803/1530 

(3.738 ± 0.008) x lO" 9 erg cnrV 

4.772 x lO" 9 erg cnrV 1 

5.196 x 10~ 9 ergcnrV 1 



3. Data reduction 

There are 1 1 .79 million events in the Chandra event file, equiva- 
lent to an average count rate of 1 17 c s -1 , in 4 ACIS-S detectors. 
The 0th order contains 3.56 million events (30%), -1/+1 orders 
contain 2.43/3.29 million events, orders -2/+2 0. 12/0. 14 million 
events, orders -3/+3 0.27/0.28 million events. We investigated 
only orders -1 and +1 since these encompassed 7.1 times more 
photons than orders -2/+2 and -3/+3 combined. 

We reduced the Chandra data with ciao version 4.4, f tools 
version 6.12 and CALDB version 4.5.0. tgextract was used 
to extract the grating order data from the event-2 file, after the 
eventfile was filtered for a particular time interval with dmcopy, 
including a workaround 5 for a bug in ciao versions 4.3 and 4.4 
concerning incorrect GTIs when time-filtering grating data with 
dmcopy. The resulting pha2 file was divided into separate or- 
ders with dmtype2split. fullgarf was used to obtain ancil- 
lary files for each time interval separately, mkgrmf for the re- 
sponse files (only one file for all time intervals). We excluded 
data around three absorption edges because the effective area 
curve there is affected by uncertainties in the model of the de- 
tector contaminant. Data were excluded in the following photon 
energy intervals: 0.52-0.55 (oxygen K edge), 0.68-0.73 (iron 
L edge) and 0.84-0.89 keV (neon K edge). Orders were com- 
bined through add_grating_orders. Despite the high count 
rates during the burst, none of the Chandra data are significantly 
piled up, due to the small frame time and large dispersion. The 
burst spectra were also independently extracted and fit using idl 
and isis, confirming all spectral fit results. 

RXTE/PCA spectra were extracted from event mode data us- 
ing ftool seextrct, filtering out artificial 'clock events'. The 
PCA spectral response was calculated with pcarsp version 
11.7.1. We find that the background model for PCU1 underes- 
timates the spectrum beyond 30 keV, where no source photons 
are expected, by about 25%. To optimize the statistical quality 
and employ as many photons as possible, we chose to exclude 
the bandpass beyond 20 keV rather than to exclude PCU1 data. 
Below 20 keV the source photon spectrum is at least four times 
brighter than the background spectrum, so that the error is at 
most 3% (at 20 keV) for the non-burst emission and less for 
the burst emission. We verified the results with those obtained 
after excluding PCU1 and do not find significantly different re- 
sults. Thus, we employ data of both active PCUs. The dead time 
fraction was calculated from standard- 1 data following the pre- 



scription at the RXTE web site 6 and was found to be up to 23.3% 
at burst peak. This fraction was taken into account when calcu- 
lating energy fluxes. All spectral bins up to 20 keV contained 
enough photons for Gaussian statistics to apply. 

4. Spectral analysis 

XSpec version 12.7.1b was employed as analysis tool. Chandra 
and RXTE burst spectra were simultaneously modeled when in- 
vestigating the continuum. RXTE data were fitted between 3 and 
20 keV and Chandra data between 0.5 and 6 keV. All spectral 
bins were grouped so that each bin contains at least 15 photons. 
No background was subtracted from the Chandra spectra, while 
the particle-induced plus cosmic diffuse background was sub- 
tracted from the RXTE spectra (employing pcabackest). The 
Chandra background is expected to be a minor part (a few c s _1 , 
see Chandra website 7 ) of the non-burst emission (117 c s _1 on 
average). We kept the normalization factor of the Chandra data 
with respect to the RXTE data to 1 . When left free, the fit proce- 
dure finds a minimum xl at a different normalization for every 
spectrum by adjusting other spectral parameters. Furthermore, 
the improvement in xl is never more than about 1%. We con- 
clude that the normalization is ill constrained. Thus, fluxes are 
calibrated against the PCA response. 

4.1. Chandra/RXTE pre-burst spectrum 

A good representative of the non-burst emission during the burst 
may be obtained from data during the 2.5 ks prior to the burst. 
The flux remains approximately constant during that time and is 
similar to the flux immediately after the burst (for a high-quality 
light curve from RXTE data, see Fig. 3). Furthermore, there is 
simultaneous Chandra and RXTE data for this interval. We de- 
rived a satisfactory empirical model for these data consisting of 
a disk black body (e.g., Mitsuda et al. 1984), a power law and 
a gaussian component, all absorbed by a single medium mod- 
eled according to the prescription in Wilms et al. (2000). The 
Gaussian component only applies to PCA data. It was also of- 
ten needed in PCA spectra of other X-ray bursts (Galloway et al. 
2008) and may be related to the collection of additional Galactic 
emission in the wide field of view of the PCA. Table 1 presents 
the best-fit model parameters. The value for /V H is within the 
range of values found previously by, for instance, Cackett et al. 
(2009) and Papitto et al. (2009). 

4.2. Burst spectrum 

We first fitted the burst data with an absorbed black body model, 
keeping the non-burst model fixed. This is the most common 
model applied to X-ray burst spectra. This model does not pro- 
vide a statistically acceptable fit to the data for the first 20 s of 
the burst, neither the combined data nor the Chandra and RXTE 
data alone. This is illustrated in Fig. 4 (top panel) which shows 
the data for 2 to 5 s after burst onset when the temperature is 
varying by no more than 10% (see below), along with the best- 
fit model; the goodness of fit is xl = 3.58 with v = 267 degrees 
of freedom. There is a strong excess at low photon energies, par- 
ticularly in the Chandra data. In fact, the Chandra data alone 
are best fitted with a power law instead of a Planck function, in 
contrast to the RXTE data. 



http : //cxc . harvard . edu 



6 http : //heasarc . gsf c . nasa . gov/ docs/xte/recipes/pca_ 
deadtime.html 

7 http : //cxc . harvard . edu/contrib/maxim/bg/index . html 
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Fig. 4. Deconvolved spectra for 2 to 5 s after burst onset. The top panel 
shows the Chandra (black) and RXTE (red) spectrum. Data points are 
indicated by crosses and the total model by the histogram. The upper 
dotted curve represents the black body component and the lower two 
curved lines comprise the model components for the pre-burst non burst 
radiation. The model is a bad fit (xl = 3.505 for v = 268), but can be 
made consistent with the data by leaving free the normalization / a of 
the non-burst model (xl = 0.963 for v = 267). In the bottom panel the 
residual spectrum is plotted after subtraction of the fitted black body 
model (crosses). The model for the non-burst spectrum is plotted as a 
histogram. For clarity Chandra data is plotted up to 3 keV and RXTE 
data beyond 3 keV (only PCU2 to extend the energy range from 20 to 
30 keV for visualization purposes). 



In order to find an explanation for the soft excess, we tested 
four spectral models. We did this in two time intervals, because 
the model needs to be applicable over all times. We chose time 
intervals 2-5 and 12-16 s, since during these times the spectrum 
does not vary much while the statistical quality is good. The first 
two models take into account inelastic scattering of photons by 
electrons in either the NS atmosphere or the accretion disk. The 
former was first put forward by London et al. (1986), further 
developed by Madej et al. (2004) and most recently calculated 
by Suleimanov et al. (2011, 2012). Hot atmospheric electrons 
harden the photons coming from below and increase the ob- 
served temperature to a value that is one to two times larger than 
the effective temperature. A soft excess remains (see Fig. 7 in 
Suleimanov et al. 201 1). The model's free parameters are the lu- 
minosity in terms of the Eddington value and the NS radius. The 
fit of this model to the first data set (see Table 2) is insufficient, 
but that may be expected because model does not formally apply 
to radiation at the Eddington limit. The fit to the second data set, 
at sub Edddington flux levels, is not acceptable either. The data 
have a broader spectrum than the model. 



Scattering by the accretion disk is calculated through the 
reflection model by Ballantyne (2004). In this model, black 
body radiation is assumed to hit the accretion disk and instantly 
photo-ionize it. The disk is modeled as a constant-density 1- 
dimensional slab. The radiation then is reprocessed by the disk 
and re-emitted into the line of sight. Above a few keV the re- 
flected spectrum is very similar to the black body spectrum. 
Below a few keV it shows a soft excess whose detail depends 
on the level of ionization. The magnitude of the soft excess 
is a strong function of the density of the disk (see Fig. 4 in 
Ballantyne 2004). The total amount of reflection depends on the 
observer's viewing angle of the accretion disk. The model's free 
parameters are the black body temperature kT, the ionization pa- 
rameter log£ (£ in units of erg cm s~') and the reflection fraction 
R. We first tested the disk model with density «h = 10 15 H-atoms 
cm -3 on the 2-5 s data, but this turns out to be unacceptable. 
A model with «h = 10 20 H-atoms cm -3 performs better and is 
able to fit the soft excess reasonably. The results are provided in 
Table 2. The reflection fraction turns out to be very large (10) for 
the 2-5 s data. It cannot be forced to smaller values by leaving 
free other parameters such as Nh- That may suggest that «h is 
even higher, although 10 23 ctrT 3 is close to the maximum value 
in standard accretion disk theory (Shakura & Sunyaev 1973). 
The strong coupling between R and «h makes it difficult to ob- 
tain a good constraint on both, though. Unfortunately, due to in- 
creasing importance of 3-body interactions, there is no model 
yet available for higher «h values. 

The third model involves the addition of a second black body 
component with a different temperature. This is based on the 
possibility that there is exists a boundary layer between the ac- 
cretion disk and the neutron star where radiation is released due 
to friction between both (Inogamov & Sunyaev 1999, 2010). Fits 
with this model are fairly good, see Table 2. We checked whether 
the small residual soft excess for the 2-5 s data is resolvable with 
a smaller value for Nh but it is not completely (xl = 1 .229 with 
y = 266for/y H = cnr 2 ). 

The fourth model is a straightforward variation of the ini- 
tial model, namely to leave free the normalization / a (with re- 
spect to the pre-burst value) of the persistent emission compo- 
nent listed in Table 1 . This simple model was recently employed 
by Worpel et al. (2013) to model successfully 332 PCA-detected 
X-ray bursts with photospheric expansion in the 3-20 keV band. 
This model performs at least as well in terms of xl as the dou- 
ble black body and reflection model, see Table 2, but on top of 
that it shows the least amount of systematic trend in the resid- 
uals for both spectra. Therefore, we chose to perform the full 
time-resolved spectral analysis with this model and note that the 
results for the primary black body component are similar as in 
the double black body model. Henceforth, we will call this the 
'fa model. 

The 'f a ' model is illustrated in the bottom panel of Fig. 4 
where the 2-5 s spectrum is shown after the fitted black body 
model is subtracted. The model of the pre-burst spectrum is also 
shown. Over the 0.5-30 keV range shown, the shape is the same 
to a fairly accurate degree. To obtain a sense of how similar the 
persistent spectrum during this interval is to that before the burst, 
we left free in addition the power law and disk black body pa- 
rameters. The best-fit power law index becomes 1 .83 + 0.05 and 
the disk black body kT m = 1.23 ± 0.64 Of 2 = 0.701 for v = 263). 
Therefore, it is somewhat harder. 

We defined time bins covering 200 s from the burst onset, 
that were sufficiently short to resolve the spectral variations, 
while maintaining adequate signal-to-noise for precise spectral 
parameters. The adopted binning, achieved via trial-and-error, 
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Table 2. Fit results when attempting four different models (columns) to account for the soft excess on 2 data sets (rows). First order LETGS spectra 
were combined and binned to 15 photons per bin. The figures show the deconvolved spectra (top panels; in units of phot cnr 2 s _1 keV -1 ) and their 
relative deviation (i.e., data divided by the predicted model value; bottom panels) with respect to the models. Black curves are for LETGS and red 
for RXTE data. The y-axes for the bottom panels have identical scales. 



I NS atmosphere model 



I Reflection model 



Double black body 



I '/ a ' model 



Time interval: 2-5 s(x'j = 3.505 with v = 268 for black body fit) 



Xi = 2.965 (y = 268) 
L/L Edd = 0.501 ±0.004 
R = 16.49 ± 0.09 km 




X; = 1.099 (v = 266) 
logf = 3.60 ± 0.03 
kT = 2.27 ± 0.03 
R = 10.0 ± 12.8 
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[ +^4^ 

h 


— Sj 




W^iMWf H 







Xt,= 1.274 (v = 266) 
kTi = 0.620 ± 0.017 
F\ =18% 

kT 2 = 2.082 ±0.013 
F 2 =82% 



Xi = 0.963 (v = 267) 
kT = 1.985 ±0.009 
F bb =63% 
/,= 16.6 ±0.6 
Ft =37% 




Xt = 1-893 (v= 148) 
L/Lem = 0.293 ± 0.002 
R = 13.76 ±0.08 km 



Time interval: 12-16 s fog = 1.792 with v - 148 for black body fit) 



XT= 1.341 (y= 147) 
kT = 1.790 ±0.007 
F bb =80% 
f a = 3.7 ± 0.3 
F fa =20% 
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Xi = 0.954 (y = 146) 
logf = 3.61 ±0.14 
kT = 1.84 ±0.02 
R = 0.38 ± 0.24 
(« H = 10 20 cm" 3 ) 




kTi = 1-593 
Fi =73% 
kT 2 = 2.82 ± 
F 2 =27% 



: 0.046 



0.34 




totaled 47 approximately logarithmically-spaced bins between 
0.1 s at burst onset and 16 s in the tail. All time bins except 
one bridge the RXTE data gaps due to the full data buffer (see 
Sect. 2). In Fig. 5 we show the results of fitting the ' f' a model. 
The fits are excellent, except perhaps between 10 and 20 s, see 
Sect. 4.3. The soft excess remains until 20-30 s after the burst 
onset. In Fig. 6 we plot / a — 1 against the burst bolometric flux. 
This shows that the increase of the persistent spectrum is present 
for burst fluxes above 10% of the peak flux and that its flux is 
not a simple constant fraction of the thermonuclear burst flux. 

The peak radius reached during the expansion interval was 
10 km, a factor of approximately 2 larger than the radius at 
touchdown. This degree of expansion is quite modest; in the 
strongest radius expansion bursts, the radius expands to at least 
10 2 times the NS radius. We note that, similar to other bursts 
from this source, the radius increases following touchdown up 
to a level comparable to the maximum reached during the ex- 
pansion. However, this increase during the burst tail is usually 
attributed to a decrease in the spectral correction factor as the 
flux decreases (e.g. Suleimanov et al. 2011) so that the photo- 
spheric radius during the expansion is significantly in excess of 
the NS radius. 

The burst reaches a peak flux of (2.50 + 0.16) x 
1Q- 7 erg cm- 2 s-' (after accounting for the increased persistent 



emission), and the bolometric burst fluence is fi, UISt = (3.41 ± 
0.04) x lO -6 erg cirr 2 over 184 s. Table 3 presents the complete 
set of derived burst parameters. One can estimate the burst recur- 
rence time according to f recU r -ax /burst /F peis , where a * 150 
(Galloway & Cumming 2006) and bolometric persistent flux 
^pers = 6 x 10~ 9 erg crrT^- 1 . This results in f recul - « 80 ks. This 
is similar to the maximum continuous data set available around 
the burst (see Fig. 2). An interesting number is the total number 
of bursts expected in the whole outburst. If we assume a Crab 
spectrum for the whole outburst, the total fluence under the esti- 
mated integral of the outburst light curve (Fig. 2) would be about 
/outburst = 2.9 x lO -3 erg cirr 2 . The total number of expected 
bursts would then be /outburst /ff/burst ~ 6. The net coverage of 
the RXTE and Chandra observations is only 2.5 d compared to 
roughly 25 d for the whole outburst. Thus, it is not unexpected 
to have detected only one burst. 

4.3. Search for absorption edges 

We extracted full-resolution orders -1 and +1 of the LETGS data 
for the first 20 s and fitted the spectrum after subtraction of the 
pre-burst spectrum with a power law and show in Fig. 7 the de- 
viations of the data with respect to the continuum model. The 
goodness of fit is xl - 0.684 (v = 678). There is no obvious 



6 



in 't Zand, Galloway, Marshall, Ballantyne et al.: A bright thermonuclear X-ray burst observed with Chandra and RXTE 



5 




0.0 E , , , 

0.1 1.0 10.0 100.0 



Time [sec] 

Fig. 5. Time-resolved spectral analysis results for the combined 
RXTE/Chandra data sets. The top panel shows the photon count rates 
in both instruments: gray for Chandra-LETGS (in units of 500 s -1 ) and 
black for RXTE-PCA (in units of 10 4 s" 1 ). 
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Fig. 6. Persistent emission flux pre-factor / a minus 1 versus bolometric 
burst flux measurements (see also Fig. 5). The gray area indicates the 
region of the Eddington limit, between that for a hydrogen and a helium 
photosphere. The dashed line shows how the trend should go if the flux 
ratio of the increased persistent spectrum to burst would be constant. 
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Fig. 7. Combined and rebinned first-order LETGS spectrum of data 
from the first 20 s of the bursts, in terms of fractional residuals with 
respect to an empirical continuum fit. 



sought-after edge visible. We determined the 90% confidence 
upper limits on the equivalent width 8 of any absorption edge by 
fixing a model for the edge at various photon energies, leaving 
free the power law parameters and edge depth while fixing Nu 
and determining the edge depths for which x 2 is 2.7 above the 
minimum value. We find the upper limit to range between 20 eV 
at 1.8-2.4 keV to 200 eV at 1.3 and 3.3 keV up to 400 eV at 5.0 
keV, see Fig. 8. 

There is a broad dip at 1.8 keV which may arise from a 
broad absorption line or edge feature. Fitting it with an absorp- 
tion line yields a centroid energy of 1 .77 + 0.04 keV, a line width 
of cr = 0.15 ± 0.05 keV and an equivalent width of -56 ± 15 eV. 
The goodness of fit is xl = 0.624 (v = 673). Fitting it with an 
absorption edge yields an edge energy of 1.29+0.01 keV and an 
optical depth of 0.29±0.06 with xl = 0.648 (v = 676). Adding 
a second edge does not really improve the fit, with^J = 0.641 
(v = 674). 



8 we define the equivalent width as the integral over energy E of 
1 - exp(-r ( ) with E the edge energy and r its optical depth 



We searched for edges in smaller time frames and found 
none. The upper limits are a factor of 4 to 5 worse at 1 s res- 
olution than at 20 s. 

In the RXTE data we analyzed the spectra with the highest 
xl values: between 6 and 8 s (this is near the touch-down point) 
and between 12 and 16 s. These spectra show spectral deviations 
which are reminiscent of the shape of an absorption edge with 
an energy between 7 and 8 keV and an optical depth of 0.2-0.3 
(equivalent width 0.7-1.1 keV). Unfortunately, these edges are 
not covered by the Chandra data so we cannot get confirmation 
from there. We repeated the time-resolved spectroscopic analy- 
sis illustrated in Fig. 5, but adding an absorption edge and find 
that the data are consistent with the continuous presence of the 
above mentioned edge, but that it is detectable only when the sta- 
tistical quality of the data are sufficient to allow detection which 
is between roughly 6 and 40 s after burst onset. 

4.4. Search for spectral lines 

We searched for absorption and emission lines in the LETGS 
data during five time intervals: 0-2 s during the burst rise, 2-5 
s during the first part of the peak, 5-9 s during the second part 
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Fig. 8. For the 20 s exposure (see Fig. 7), this graph shows the upper 
limits of an absorption edge as a function of photon energy. 



Wavelength (Angstrom) 

Fig. 9. Chandra-LETGS burst spectrum between and 7 s after burst 
onset, in wavelength domain. The histogram represents the data and the 
dotted line a continuum fit. Top panel is the spectrum, bottom panel the 
deviation with respect to the fitted model in units of a. Note the feature 
at 11.9 A. 



of the peak, 9-12 s during the first part of the tail and 12-16 s 
during the second part of the tail. A 4. Oct detection threshold 
was employed. This is equivalent to one chance detection in all 
20,000 trials. The number of trials is determined as follows: 200 
independent tests per 2-12 A spectrum with a LETGS resolution 
and binning of 0.05 A; a factor of 2 for five successive 2 times 
binned-up versions of the spectra; and a factor of 5 for the five 
time intervals mentioned above. In order to have a 10% chance 
of a false positive (in either absorption or emission), a one-tailed 
test would require a 1 in 20,000 chance for detection in random 
data, which is about 4.0 sigma. We also searched in 4 different 
time intervals (0-7, 2-4, 4-6, and 7-30.5 s) that overlap the above 
time intervals. 

The strongest candidate line that we find is an emission fea- 
ture that we find in the 0-7 s time interval at 1 1 .9 A (or 1 .04 keV), 
see Fig. 9. This feature is 4.7cr above the continuum in a binning 
of 0.11 A. Given that we searched in additional time intervals, 
the significance should be diminished to 4.6<x. The equivalent 
width is about 4 A. 

5. Timing analysis 

We performed a timing analysis of the RXTE data of the burst, 
by generating Fourier power spectra for 10-s data stretches every 
1 s, employing the event mode data with 125 yus resolution. We 
limited the analysis to channels 5-37 (2-10 keV), where the ra- 
tio of the burst photon rate to the background rate is highest. 
We detect burst oscillations with a burst fractional amplitude 
of 2-5% rms in the tail. This value is calculated without tak- 
ing into account the contribution of the accretion powered pul- 
sations that might still be present at the time of the burst. This 
correction is relatively small since the accretion powered pul- 
sations have pulsed fractions of 3-4% rms prior and after the 
occurrence of the burst (see for example Patruno et al. 2012). 
Furthermore the photon count rate detected when burst oscil- 
lations are seen is between two and four times the count rate 
received prior to the occurrence of the burst, meaning that the 
oscillations seen are truly of nuclear origin (see e.g., Eq. 11 in 



Watts et al. 2005). Similar considerations apply to the burst oscil- 
lation phase, which are negligibly influenced by the presence of 
accretion powered pulses. Burst oscillations were marginally de- 
tected (between 2 and 3 sigma for 1 s long trains) at 400-401 Hz 
during the burst rise with fractional rms amplitudes of up to 5%. 
The formal 90%-confidence upper limit on the fractional rms 
amplitude is 8%. This is consistent with measurements of pre- 
vious bursts from SAX J1808.4-3658 with RXTE (Chakrabarty 
et al. 2003; Bhattacharyya & Strohmayer 2006) where ampli- 
tudes were seen between 5 and 25% in 0.25 s long trains. Our 
data is less sensitive because of the smaller number of active 
PCUs (2 instead of 4). 

The time resolution of the ACIS-S CC mode data is 2.85 
ms (351 Hz). This implies that the burst oscillation signal (401 
Hz) cannot be easily resolved in the Chandra data. Simulations 
show that to detect a sinusoidal pulsation at 401 Hz in a 1 s 
stretch of data with a resolution of 2.85 ms requires a relative 
rms amplitude of at least 25% to be significantly detected as an 
alias at 401-351=50 Hz which is an order magnitude larger than 
ever detected for SAX J1808.4-3658. 



6. Discussion 

The 2011 outburst of SAX J1808.4-3658 is the seventh since 
its discovery in 1996 (in 't Zand et al. 1998), and comes 3 years 
after the previous outburst. Galloway (2008) found that the onset 
time of the earlier outbursts followed a quadratic relation with 
time, which predicted the time of the previous outburst to within 
lid. However, the 2011 outburst occurred 197 d earlier than 
predicted, suggesting that the long-term mass transfer rate may 
be increasing. We estimated the total fluence for the outburst as 
=s 2.9 x 10~ 3 erg cm 2 , which is significantly lower than for the 
first four outbursts (Galloway 2006). Thus, it seems likely that 
the long-term accretion rate remains substantially below that of 
the first decade of observations. 
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Table 3. Burst parameter values, extracted from RXTE data alone for 
easy comparison with other bursts (e.g., Galloway et al. 2008). Times 
are non-barycentered. 



Parameter Value 

Onset time MJD 55873.91320 (barycentered) 

2011 Nov 8 21:55:00.50 UTC 
563407175.90 s RXTE MET 
437176836.46 s Chandra MET 

Rise times 2.5 ± 0. 1 s (25 to 90% of peak) 

6.3 ±0.1 s(0 to 100%) 
0.1 ±0.1 s(0to 70%) 

e-folding decay time 12.0 s 

Time scale (fluence/peak flux) 12.9s 

Duration «500 s (down to 0.03% of peak flux) 

Unabs. bol. peak flux (2.56 ± 0.15) x 10~ 7 erg cnrV 1 

Unabs. bol. fiuence (3.41 ± 0.04) x 10~ 6 erg cnr 2 

3-20 keV persistent flux 1.72 x 10~ 9 erg cirrV 

0.5-60 keV unabsorbed flux 6.1 x 10~ 9 erg cirrV 1 



6.1. Peak flux and burst oscillations 

The thermonuclear burst detected from SAX J1808.4-3658 on 
2011 November 8 is the brightest ever detected with Chandra, 
and the second brightest with RXTE. Within the 1907 ks expo- 
sure time on SAX J1808.4-3658 during the complete RXTE 
mission, this is the ninth burst (e.g., Galloway et al. 2008). 
The second brightest burst has a peak flux of (2.40 + 0.03) x 
10~ 7 erg cm s~' (this is burst no. 4 in Galloway et al. 2008, 
note that this publication lists peak fluxes before dead time cor- 
rection) that is only slightly fainter than the brightest. Our burst 
is otherwise fairly ordinary for SAX J1808. 4-3658, including its 
timing behavior. Galloway & Cumming (2006) find that bursts 
from SAX J1808.4-3658 are due to flashes in a hydrogen- 
depleted layer and burn primarily helium. 

The Eddington limit of a solar-composition photosphere on 
the surface of a 10 km radius NS of mass 1.4 M at a dis- 
tance of 3.5 kpc corresponds to a flux of 1.4xl0~ 7 erg cm _2 s _1 , 
for a helium-rich photosphere 2.6x10 7 erg cm ^s 1 . The peak 
flux for the burst from SAX J1808.4-3658 is within this range, 
clearly indicating that the energy production rate due to the ther- 
monuclear flash is consistent with the Eddington rate. It is per- 
haps worth noting here that the distance estimate of 3.5 kpc 
for SAX J1808.4-3658 comes from a comparison of the mea- 
sured burst energetics and recurrence times with the predictions 
of theoretical ignition models (Galloway & Cumming 2006), not 
merely from the peak flux of previous PRE bursts, as is the case 
for the majority of LMXB distances. While the time dependence 
of the black body temperature and normalization is also charac- 
teristic for an Eddington-limited burst, with the typical ' touch- 
down' peak at 10 s (see Fig. 5), the radius expansion is limited 
to at most a factor of about 2. This is typical for moderate ex- 
pansion (in 't Zand & Weinberg 2010). There is no evidence for 
the expulsion of a shell, such as a (short) precursor. Ergo, this is 
not a superexpansion burst (see Sect. 1). 

The accretion flux prior to the burst translates to a 0.5-20 
keV luminosity of 4.5xl0 36 erg s _1 for a distance of 3.5 kpc, or 
64 times smaller than the burst peak in the same bandpass. 

The burst oscillations detected in this burst are very similar 
to those detected in the RXTE data of previous bursts. The am- 
plitude compares well with the fractional rms amplitude of 2-4% 
in previous bursts (Watts 2012; Patruno & Watts 2012). 



6.2. Broad-band spectrum 

The Chandra data are interesting for the peculiar behavior of the 
burst spectral continuum. As illustrated in Fig. 4, data below 3-4 
keV draw a rather different picture of the burst spectrum than 
above. The sub-3 keV data reveal a strong soft excess above the 
black body. This is not the first time that a soft excess is observed 
in a burst spectrum, although it seems to be a rare phenomenon, 
possibly because many bursters have too much interstellar ab- 
sorption to allow detection of a soft excess. We count about 250 
X-ray bursts that have been detected with sub 3 keV coverage by 
XMM-Newton, Chandra and Swift-XRT (e.g., Kong et al. 2007; 
in 't Zand et al. 2009; Galloway et al. 2010). Most of these do not 
show photospheric expansion. Only Asai & Dotani (2006) and 
Boirin et al. (2007) report on strong soft excesses found in some 
bursts from the eclipsing LMXB EXO 0748-676. They propose 
that this is due to a decrease in absorption by the circumstellar 
medium arising from the photo-ionization of that medium by the 
burst photons, thereby allowing more accretion disk flux to reach 
the X-ray detector. This explanation may be plausible for a high- 
inclination system (i = 75°-83° for EXO 0748-676; Parmar 
et al. (1986)), but for SAX J 1808.4-3658 the inclination angle is 
suggested to be low (< 67°; Deloye et al. 2008), so there is little 
local obscuring material affecting the spectrum. We note that the 
two published Chandra HETGS grating observations of 33 X- 
ray bursts (including 4 with photospheric expansion; Thompson 
et al. 2005; Galloway et al. 2010) could have found soft excesses 
of similar magnitude as in SAX J1808.4-3658, had they been 
present. 

The net burst spectrum is not consistent with the NS atmo- 
sphere model of Suleimanov et al. (201 1). There is better agree- 
ment with the model for reflection against the accretion disk and 
even better with the simple double black body model. Mildly 
better overall is the agreement with the 'f a ' model, because 
this models appears to have the smallest systematic trends in its 
residuals (see bottom panels of figures in Table 2). In a study 
of 332 radius-expansion bursts seen with RXTE, Worpel et al. 
(2013) modeled the 3-20 keV spectra in the same '/„' manner. 
They found as a general feature that f a is significantly in excess 
of 1, with the highest values arising from the bursts observed at 
the lowest pre-burst fluxes. Our combined Chandra-RXTE anal- 
ysis of one burst from SAX J1808. 4-3658 is consistent with 
those results and shows additionally that this behavior extends 
to lower energies, giving rise to a large flux excess below 3 ke V 
for a system with low Nh (see Fig. 4). However, we note that the 
difference in performance with other models is not much better. 
It is in the spirit of Occam's Razor (with only 1 additional free 
parameter) that we provide the full time-resolved spectroscopic 
analysis in terms of the simple '/„' model, but other models are 
almost equally justified. 

Worpel et al. (2013) suggest that the increase of the non-burst 
emission is caused by the Poynting-Robertson (PR) effect, which 
removes angular momentum from the inner accretion disk due to 
radiation drag, thereby increasing the accretion rate during the 
burst. If the accretion rate in SAX J1808.4-3658 increased by 
a factor of 20 during the burst, based on our maximum best-fit 
value of f a , it would come close to 50% of the Eddington limit. 
Significant changes in the persistent spectral state might be ex- 
pected, analogous to the so-called spectral state changes of the 
persistent emission that are attributed to changes in the geome- 
try of the accretion disk (e.g., Done et al. 2007, and references 
therein). Since the persistent spectrum of SAX J1808.4-3658 at 
the time of the burst is quite hard, a softening of the spectrum 
seems most likely. However, it is not obvious that a typical hard- 
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to-soft spectral state change can take place in response to an in- 
crease in accretion rate sustained for < 20 s, or even if the effects 
of PR drag would produce the same observable effect. Our data 
suggests that, if there is a change, the persistent spectrum instead 
becomes harder during the burst. 

We note also that if the accretion rate increases by a fac- 
tor of 20 over 10 s, the additional matter must be that which 
would otherwise have been accreted over the next 200 s, and one 
would thus expect a suppression of the accretion rate following 
the burst. This would be the hallmark of the PR effect taking 
place. We found no evidence of f a values smaller than 1, out to 
500 s after the burst. However, the recuperation time of the disk 
may be longer if the disk viscosity is relatively small (Walker 
1992), so that the suppressed accretion rate would be spread out 
over a longer interval, and hence undetectable. 

The pre-burst spectrum likely includes a component from 
inverse Compton scattering of undetectable soft photons (with 
energies < 0.5 keV for this burst) by an accretion disk corona 
(ADC; White & Holt 1982; Begelman et al. 1983; Begelman 
& McKee 1983). An increase in the persistent spectrum would 
therefore require either a larger ADC, or an increase of soft seed 
photons into the ADC. The burst cannot produce a larger ADC 
through photoionization, as the temperature of the plasma can 
only be raised to its Compton temperature which is « kT for 
a blackbody spectrum. This temperature is much lower than the 
> 10 keV temperature required to explain the high-energy extent 
of the persistent emission. Likewise, since the local dissipation 
in the accretion disk scales as R~ 3 (where R is the radius along 
the disk), expansion of the ADC outwards through a wind or 
outflow would not be able to explain the significant increase in 
the ADC emission. This constraint could be circumvented if the 
disk and ADC were expanded inwards by the burst, as through 
Poynting-Roberston drag (see discussion above). 

Perhaps the simplest explanation for an increase in the per- 
sistent spectrum during the burst is a larger influx of seed pho- 
tons produced by reprocessing in the accretion disk. As dis- 
cussed in Sect. 4.2, high density reflection models were able to 
explain a large fraction of the strong soft excess seen early on in 
the burst. This emission, which extends below 0.5 keV, is dom- 
inated by thermal bremsstrahlung produced by the X-ray heated 
outer layers of the disk (Ballantyne 2004; Ballantyne & Everett 
2005), and is a natural source of additional seed photons to the 
corona. This explanation has several advantages: first, the ab- 
sorption and re-emission of the burst will all occur on not much 
more than twice the light-crossing time of the system, and would 
explain why the enhanced persistent emission is seen even in the 
earliest phases of the burst. Similarly, Figs. 5 and 6 illustrate that 
the increase in the persistent flux falls rapidly as the illuminat- 
ing blackbody fades, which would be a natural consequence of 
less reprocessing in the inner disk. Finally, if reprocessing ac- 
counts for most of the soft excess, then this reduces the values 
of f a required to fit the spectra. Nevertheless, this idea remains 
just informed speculation, and further study of the broadband 
spectra of Type 1 X-ray bursts, as well as a concentrated theoret- 
ical effort into understanding the interaction of a burst with the 
surrounding accretion flow, is needed to understand the physical 
processes at work in these environments. 

6.3. Evidence for spectral features 

The high spectral resolution data provided by the 
Chandm/LETGS show no evidence for the sought-after 
absorption edge with upper limits in the equivalent width of at 
best 20 eV. The residuals between 1 and 2 keV with respect to 
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Fig. 10. For the 20 s exposure (see Figs. 7 and 8), this graph shows 
the upper limits of an absorption edge at 2.67 keV (lower horizontal 
dash-dotted line) and 3.48 keV (upper horizontal line), compared to the 
predictions for the K-edge of hydrogen-like 28 Si (dotted curve), 32 Ca 
(solid curve) and 56 Fe (dashed curve), for abundances as determined 
from the HeO.l model (pure helium ignition at 0.1 Mem) in Weinberg 
et al. (2006). The range where the upper limits apply are set by the 
range of black body color temperatures measured (2-3 keV) multiplied 
with a maximum possible color correction of 0.6 to arrive at effective 
temperatures (e.g., Suleimanov et al. 201 1). The upper limits would be 
3 to 4 times better if the edge energy would be gravitationally redshifted 
by a factor of 1.31. They are 4.5 times worse for an exposure time of 1 
s. 



a continuum fit (Fig. 7) are not clearly identifiable with a single 
absorption edge or two. A single edge would have an energy 
of 1.29 keV. For a gravitational redshift anywhere between 
1 + z = 1 for far away from the NS and 1 + z = 1 .3 1 on the 
surface (for a NS with mass 1.4 M Q and radius 10 km), the 
expected range of edge energies is 1.3-1.8 keV. This is not 
identifiable with a K or L-edge of any likely element. The 
redshifted K-edge of hydrogen-like Si comes closest at 2.0 ke V. 
Only K-edge energies of Mg and Al and L-edge energies of 
Ga, Ge, As, Se and Br would fit. Also, the residuals in this 
energy range are likely to be affected by a sharp change in the 
effective area arising from the K-edge of the neutral silicon 
CCD material. 

For a 20 s exposure, the upper limit on the equivalent width 
ranges between 42 eV at the Si K-edge (2.67 keV) and 124 eV at 
the Ca K-edge (3.48 keV). These are the two most abundant ele- 
ments expected for pure helium ignition according to Weinberg 
et al. (2006) with the largest predicted absorption edges. These 
upper limits are compared with predicted values in Fig. 10. The 
predictions are for ignition in a pure helium layer like for S AX 
J1808.4-3658 (Galloway & Cumming 2006). The predictions 
are at (for Si) or below the measured upper limits. Thus, these 
measurements are consistent with predictions. We note that the 
tentative edge detected in the RXTE data at about 7.9 keV has an 
equivalent width of 0.7-1.1 keV and is, thus, much higher than 
the prediction for the Fe-K edge. On the other hand, the picture 
that this edge draws is not consistent with the expectation. The 
edge energy and optical depth are consistent with being con- 
stant while the temperature of the photo-ionizing radiation field 
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changes between 3 and 1 keV. This is different from the features 
seen by in 't Zand & Weinberg (2010) and suggests one should 
be careful in associating this feature with matter close to the NS. 

We find marginal evidence for an emission line at 11.9 A 
with an equivalent width of 4 A. This wavelength resides in the 
range of rest wavelengths of a number of Ne lines and is slightly 
blue-shifted from the Fe L line complex. We believe the line is 
too marginal to draw any further physical conclusions. However, 
it is interesting to note that recently a strong emission line has 
been detected (Degenaar et al. 2012) at about the same energy 
in a long superexpansion burst, with an equivalent width that is 
about 3 times smaller than the feature we see, but for a much 
longer duration (10 min). 

Although this burst is Eddington limited, it does not ex- 
hibit strong photospheric expansion. Such expansion will probe 
the predictions at lower black body temperatures in Fig. 10. 
Thus, the sensitivity would be better. We suggest to focus fu- 
ture searches for absorption edges on burst with those, ergo the 
'superexpansion' bursts, since these have not yet been detected 
at high spectral resolution. The most opportune sources to search 
in are 4U 1820-30 (in globular cluster NGC 6624) and 4U 1724- 
30 (in globular cluster Terzan 2), because they show bursts every 
few hours to days (4U 1 820-30 in the low/hard state) with at least 
25% of the time superexpansion (in 't Zand & Weinberg 2010; 
in 't Zand et al. 2012). 

7. Conclusion 

A simultaneous Chandra-LETGS/RXTE detection of a very 
bright thermonuclear X-ray burst with moderate photospheric 
expansion and low Nu presents a unique opportunity to study 
the 0.5-30 keV spectrum of such an event, with large sensitivity 
and high spectral resolution below 6 keV. We find that 

• the spectrum show a strong deviation from black body radia- 
tion, particularly at low energies where it exhibits a soft ex- 
cess that exceeds the black body by an order of magnitude. 
The non-Planckian component has a similar spectral shape as 
the pre-burst spectrum that is powered by accretion, suggest- 
ing that the accretion flow changes during the burst, probably 
due to the near-Eddington flux from the burst. Exactly how 
that happens is difficult to determine, mostly because there is 
no detailed model yet available for the dynamic and radiative 
interaction of a burst with the pre-existing accretion flow, but 
most likely it is related to an increase of seed photons in the 
ADC. Future work should include a hydrodynamic model of 
the interaction between X-ray bursts and accretion disks and 
a comprehensive of study of f. d values over bursts with and 
without photospheric expansion; 

• no unambiguous spectral features are detected. The upper 
limits for absorption edges are, however, consistent with 
the theoretically predicted equivalent widths due to ejected 
ashes. It is expected that edges are deeper for bursts with 
photospheric superexpansion. Therefore, future searches with 
Chandra should preferably concentrate on those. This is also 
underlined by the recent finding by Degenaar et al. (2012). 
Regarding absorption edges in bursts with moderate pho- 
tospheric expansion, investigating superbursts with current 
grating spectrometers would be interesting and, in the farther 
future, the proposed LOFT mission with a 10 m 2 silicon drift 
detector (Feroci et al. 2012) has the right combination of sen- 
sitivity and spectral resolution to make progress. 
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